Harmonic oscillator model for the atom-surface Casimir-Polder interaction energy 
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In this paper we consider a quantum harmonic oscillator interacting with the electromagnetic 
radiation field in the presence of a boundary condition preserving the continuous spectrum of the 
field, such as an infinite perfectly conducting plate. Using an appropriate Bogoliubov-type trans- 
formation we can diagonalize exactly the Ifamiltonian of our system in the continuum limit and 
obtain non-perturbative expressions for its ground-state energy. From the expressions found, the 
atom-wall Casimir-Polder interaction energy can be obtained, and well-know lowest-order results 
are recovered as a limiting case. Use and advantage of this method for dealing with other systems 
where perturbation theory cannot be used is also discussed. 

PACS numbers: 12.20.Ds, 42.50.Ct 
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I. INTRODUCTION 

Many physical phenomena, including Casimir and 
Casimir-Polder forces or the Lamb shift, are related to 
the energy shift of the ground state of a system with 
a discrete spectrum interacting with a field with a con- 
tinuum spectrum [1]. Casimir and Casimir-Polder in- 
teractions have recently received much attention in the 
literature being a direct manifestation of the quantum 
properties of the electromagnetic field, in particular of 
the existence of vacuum fluctuations, even at a macro- 
scopic level [2]. They also have raised great interest for 
their importance in applications to nanotechnological de- 
vices such as micro-electromechanical systems (MEMS) 
[3]. 

The Casimir-Polder atom-surface potential is a long- 
range interaction between a neutral atom or molecule 
(or in general a polarizable body) and a conducting wall 
in vacuo [4], and it is usually considered for systems in 
their ground state. This potential arises from the inter- 
action of the atom with the vacuum fluctuations of the 
quantum electromagnetic field modified by the presence 
of the conducting wall: this yields an interaction energy 
depending on the atom- wall distance and thus a force on 
the atom. This force has a pure quantum origin, and for 
ground state atoms can be obtained by a second-order 
perturbative calculation [5]. 

A similar force exists also for atoms in an excited 
state, with a different dependence on the atom-wall dis- 
tance due to the presence of atom-field resonances [6, 7]. 
Perturbative calculations of Casimir forces for excited 
states may encounter difficulties due to resonance di- 
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vergences. There are indeed some controversies in the 
literature when Casimir-Polder forces are considered for 
systems in their excited state, in particular when a res- 
onant energy transfer between two identical atoms can 
occur [8, 9]. Also perturbative calculations of atomic 
radiative energy shifts encounter difficulties related to 
vanishing energy denominators [10]. Time dependence 
of the atom- wall Casimir-Polder interaction leading to a 
dynamical Casimir-Polder effect has been also recently 
considered in the literature [11, 12]. 

Nonperturbative calculations may avoid the difficul- 
ties arising for diverging energy denominators and, at 
the same time, be useful for obtaining Casimir-like in- 
teractions for other physical systems where the coupling 
with the field could be large or when the density of pho- 
ton states is large. For example, a large photon den- 
sity of states occurs near the band-edge of a photonic 
crystal [13], and this may require nonperturbative ap- 
proaches. In condensed matter physics, the density of 
states of one-dimensional systems such as a quantum wire 
has the well-known Van Hove singularity near the edges 
of the band, and its role in enhancing decay rates [14] 
or charge transfer from one or multiple impurities to the 
wire [15, 16] have been recently investigated with non- 
perturbative methods. An analogous effect investigated 
in the last years is the increase of the Lamb shift of an 
excited level of an atom in a photonic crystal [17]. Non- 
perturbative methods with a dressed-state approach have 
been used for describing the radiative spontaneous decay 
from an excited state [18, 19] or the evolution of an un- 
stable system [20]. 

In a previous paper we have used a method based on 
Bogoliubov-like transformations in order to obtain the 
Casimir-Polder interaction between two atoms, where the 
atoms were modeled as three-dimensional harmonic os- 
cillators [21]. This method based on Bogoliubov-type 
transformation was originally introduced in Ref. [22] for 
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a one-dimensional scalar field and applied to study the 
dynamics of an oscillator interacting with a scalar field 
[23]. In this paper we shall develop a similar method 
for treating the atom- wall Casimir-Polder interaction ap- 
proximating the atom as a quantum harmonic oscillator 
and diagonalizing the interacting Hamiltonian of the sys- 
tem. We obtain two exact expressions (within dipole ap- 
proximation) for the ground-state energy, from which the 
oscillator-wall Casimir-Polder energy can be obtained. 
One of the two exact energy shifts obtained is in terms 
of the resolvent, and the other in terms of the function 
resulting from the factorization problem of the resolvent. 
The known second-order expression of the force is re- 
covered as a limiting case from both expressions. Also, 
we stress that the validity of our method is not limited 
to this case and can be applied to any boundary condi- 
tion that yields a continuous spectrum of the field. In 
fact, the boundary condition enters only in the specific 
form of the mode functions of the field, and consequently 
in oscillator-field coupling constant: our diagonalization 
procedure of the Hamiltonian is valid for a general ex- 
pression of the coupling constant. We specify its form 
only at the very end of the calculation. Moreover, hav- 
ing obtained a diagonal form of the Hamiltonian, could 
make easier dealing with different states of the interact- 
ing system, such as thermal states, excited states or in 
general non-equilibrium states. 

Examples of possible applications of our method 
to other physical systems, such as impurities in low- 
dimensional systems in condensed matter physics or dis- 
persion interactions between atoms or quantum dots in 
a photonic crystal, are also briefly discussed. 

The outline of this paper is as follows. In Section 
II we introduce our Hamiltonian model and the im- 
proper Bogoliubov-type transformation that diagonalizes 
the Hamiltonian of the oscillator interacting with the ra- 
diation held inclusive of the boundary condition; we then 
obtain two alternative (but equivalent) exact expressions 
for the ground-state energy shift. In Section HI we apply 
the results of the previous Section in order to obtain the 
oscillator-surface Casimir-Polder potential energy, recov- 
ering the known lowest-order result as a limiting case; we 
also show the importance of the resonance poles in the 



expressions obtained for the energy shift of the ground 
state. In Section IV we summarize our results and dis- 
cuss their possible relevance in different physical systems 
where strong coupling constants or high density of states 
require the use of nonperturbative methods. 



II. THE HAMILTONIAN MODEL AND THE 
BOGOLIUBOV TRANSFORMATION 

We consider a quantum harmonic oscillator with fre- 
quency ujQ = cko, described by annihilation and creation 
operators a and a^, interacting with the quantum electro- 
magnetic field described by the annihilation and creation 
operators 6kj and fc]^^-. The presence of the boundary 
conditions is mathematically included in the form of the 
coupling constant /kj and of the field mode functions 

fk,(r). 

The Hamiltonian of our system is 
H = hckoa^a + hckb]^-buij 

ki 

+ ^hj{rA){a + a^){buj-bl) , (1) 

kj 

where the third term is the atom-field interaction in the 
multipolar coupling scheme within dipole approximation. 
Phases have been chosen such that the coupling constants 
/kj, which are evaluated at the atom/oscillator's position 
r^, be pure imaginary; they are given by 

/kj (r ) = « Y M • fkj (r) , (2) 

where {x is the matrix element (real) of the oscillator 
dipole moment and V the quantization volume. 

In (2), fkj(r) arc the mode functions of the electromag- 
netic field and take into account possible boundaries. If 
a conducting plate at ^ = is present, they can be ob- 
tained from those of a perfectly conducting cubical cav- 
ity of volume V = with walls {-L/2 < x, y < L/2, 
0<z<L) 



{fkj)x = \/8(ekj)xCos 
(fkj)y = V8{ekj)ysm 
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cos(kzz) , 



(3) 



where k^ = Itt/L, ky = tott/L, fc^ = wk / L (l^m^n = order to switch from the cavity to the wall at z = 0, 
0, 1, 2, . . . ) and ekj are polarization unit vectors [25]. In at the end of the calculations one has to take the limit 



3 



L — > 00. 

Our Hamiltonian (1) also includes counter-rotating 
terms yielding virtual processes. It is well known that 
they are essential for correctly obtaining radiative pro- 
cesses such as the Lamb shift or Casimir forces for 
ground-state systems [26]. This is the great advantage 
of modeling the atom or polarizable body as a harmonic 
oscillator: it allows an exact diagonalization of the total 
Hamiltonian, even maintaining counter-rotating terms. 
Although at the end we shall focus our discussion on the 
specific case of the atom-surface Casimir-Polder force, 
our method is more general and can be applied to any 
boundary conditions, provided the field spectrum is con- 
tinuous. In fact, the specific boundary condition enters 
only in the form of /kj in (1), and we shall specify it only 
at the very end of our calculations. 

Following a similar method to that used in [21, 23], we 
introduce new bosonic operators related to the old ones 
by the following Bogoliubov-like transformations 



Ski = tija + r*^,a^ + Y,{Tt/*b^^,+Rl^'*bl,,), 
k'i' 

Blj = ik.at + r^^a + ^ (r^/'bl^, + B^fb^,,,) . (4) 



k'j' 



We determine the coefficients in (4) by requiring that 
the transformed operators satisfy free-field commutation 
relations with the Hamiltonian 



[iJ, Bkj] = —hckBi^j 



(5) 



and its Hermitian conjugate relation. In (5) the frequen- 
cies associated to the new operators are the same as the 
old ones. This is strictly true in the continuum limit 

{V = —i' oo), and in the discrete case there arc correc- 
tions 0(1 /V). Even if for simplicity we keep the notation 
with discrete sums, it is understood that we shall take the 
continuum limit at the end of our calculations. 

Transformations (4) are indeed improper Bogoliubov 
transformations, because there is not a one-to-one corre- 
spondence between new and old operators: in the; diag- 
onalized form, the discrete state disappears in the con- 
tinuum [22]. This is similar to what happens also in the 
Friedrichs model, describing a two-level system coupled 
with a field continuum in the rotating wave approxima- 
tion [27, 28]. Use of (4) in (5) yields a set of coupled 
equations for the coefficients, from which we obtain 



'5kk'<J,y + ^^^ 
k ~ ko 



k + ko 



*k7 



k' + k 



(6) 

7G(fc)/k,/k',. , (7) 
(8) 

(9) 



with the resolvent given by 

{Gik))-' = kl-l 



1 ^ ikpk' 2 
2^ k"^ - A;'2 •'^'1' 



{he) 



k'j' 



If the coupling constant is such that 



(hcy^ k' 

k'j' 



(10) 



(11) 



the resolvent G{z) is analytic in the complex z plane 
apart for a cut on the real axis. The resolvent can be an- 
alytically continued on the second Riemann sheet, where 
it can have complex poles related to energy shifts and 
decay rates [23, 24]. 

After some algebraic calculations, we can also obtain 
the inverse relations expressing the old operators in terms 
of the new ones 



at = ^(A'^^B+.+r'^^Bki) , 

kj 



k'i' 



where the coefficients are given by 



^kj 



k'i' 

^k/ 



k'i' 
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(fc + fco)G(fc)-/ki, 
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-Gik) 



{he) 



(12) 
(13) 

(14) 



2/ki/k'j' , (15) 



^ki ^ _Q,_k,)G*{k)-fl^. 
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he 



(16) 
(17) 



Substitution of (12), (13) and their Hermitian conju- 
gate with the coefficients above into the Hamiltonian (1), 
after lengthy algebraic calculations, allows to obtain H 
in a diagonal form in terms of the new operators plus an 
energy shift Eg as 



H = ^hekBl.Bkj + Eg . 

ki 

The energy shift Eg is given by 



(18) 
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he 
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JkTW 



(13) 



This expression for the ground-state energy shift is ex- 
act, and valid at any order in the coupling constant. As 
we will see in the next section, when using this expression 
much attention must be paid because of the singularities 
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of the resolvent G{z) in the complex plane. For exam- 
ple, the second term on the RHS of (19), in the limit of 
vanishing coupling constant has two second-order poles 
at A; = ±fco; these second-order poles separate each other 
when higher-order terms in the coupling are included. 

Just for the purpose to obtain Eg , there is another way 
similar to that used in [23], main difference in our case 
being the presence of the boundary condition and having 
a vector rather than a scalar field. From the general 
theory of Bogoliubov transformation, the dressed ground 
state |0)d of the full Hamiltonian H can be obtained by 
a transformation on the bare ground state (ground state 
of the unperturbed Hamiltonian Hq) as 



(20) 



where |0) is the non- interacting ground state and the 
operator y is a quadratic form of the creation operators. 
In our case, we have 



V = 



6^ 6t 



1 ^ kofkjTjjk) + 



2{ncko + Eg) 



(21) 



where ri{z) is the solution of the factorization problem 
for the resolvent [22] 



G(z) = v{zM-z) 



(22) 



The fimction ri{z) in (22) is analytic in the complex z 
plane apart from a cut in the negative real axis. 
The average energy is then obtained by 



Eg= d<0\H\0>d=<0\e 



(23) 



After some algebraic calculations, using (21) we can ob- 
tain the expression for the shift of the ground-state en- 
ergy shift 



Eg = 



ko 



1 



ko 



E„ he ^ 



hi f V{k) 



(24) 



he ^9 



This is an exact equation for the ground-state energy 
shift Eg. We are interested in solutions of (24) that are 
analytic in the coupling constant fkj, that is such that 
Eg-^0 for fuj 0. 

We have so obtained two different expressions for the 
energy shift of the ground state, given by (19) and (24): 
the former is in terms of the resolvent G{z) and the latter 
in terms of the function 77(2) resulting from the factor- 
ization problem of the resolvent. In the next Section 
we shall use these expressions to recover, at the lowest 
(sccond)-order limit, the Casimir-Polder energy between 
an atom/oscillator and a conducting plate. 



III. THE GROUND-STATE ENERGY AND THE 
ATOM-SURFACE CASIMIR-POLDER 
POTENTIAL ENERGY 



We will now show that the expressions of the ground- 
state energy (19) and (24) can be used to obtain the 
Casimir-Polder interaction between an atom, modeled 
as a harmonic oscillator, and a surface (or any bound- 
ary condition preserving the continuous character of the 
modes of the electromagnetic field) . We need to use the 
appropriate coupling constant /kj. For a perfectly con- 
ducting infinite wall, the coupling constant is given by 
(2) with the normal modes in (3). The energy depends 
on the distance between the oscillator/atom and the con- 
ducting wall due to the presence of the functions /kj(rA) 
evaluated at the atomic position: minus its derivative 
with respect to ) ) yields the oscillator-wall Casimir- 
Polder force. Up to this point this result is exact within 
the dipole approximation. 

We shall now show that the lowest order approximation 
of Eg will correctly reproduces the well-known second- 
order atom-wall Casimir-Polder interaction energy. As 
mentioned in the previous Section, approximation of (19) 
must be done carefully due to poles of the resolvent G{z) 
in the second Riemann sheet. We will first show that 
the term in the second line, which apparently gives a 
contribution starting from fourth order only, indeed con- 
tributes to the second-order energy shift due to the res- 
onance poles of G{z). In fact, when the sum over k' is 
performed in the continuum limit, the poles of G{z) in 
the second Riemann sheet yield a factor in the denomi- 
nator proportional to the square of the coupling constant 
(more specifically to the decay rate of the oscillator) , and 
this makes this contribution a second-order one. We will 
now show this in the lowest-order approximation. 

We have the following relation in the continuum limit 
between G+{k) = G{k -\- iO) and G-{k) = G{k - iO), i.e. 
the extensions of G{k) in the upper and lower complex 
plane, respectively, 



(G+(fc))-i-(G-(fc))-i 
1 V 



-Awiko 



[hcf (2 



7r)3^/ 



1 /kj \^ k"^ ■ (25) 



Using this relation in the sum/integral over (k'j') in (19), 
we obtain 



k'j' 



k + 2k' 



{k + k'Y 



1 2fco 
4fc2 (fc + kof 
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where we have used (at the lowest order) 



G-{k')-G+{k') 
1 / 1 



1 



ko + k' \ko- k' + iO ko - k' - iO 

-2wi-^5{ko - k') . 
ko + k' ^ 



(27) 



This explicitly shows that (27) is indeed of zeroth-order 
in the interaction due to the resonance poles of the resol- 
vent. Substitution of (26) into (19) gives then the energy 
shift approximated to the second order (continuum limit 
is understood) 



-y 



he 



_ 27r 



/kj(l-A) 

k + ko 

k 



V ^k + ko 

kj 



M-fkj(rA) 1^ 



(28) 



that coincides with the results that can be obtained by 
perturbation theory [6, 7]. 

The same approximated result can be also obtained 
from the expression (24) of Eg. By substituting Eg in 
the RHS of (24) with its zeroth order value, that is zero, 
we obtain 



4^^ = ^E^^°Hfc)/^.- 



(29) 



kj 



where r]'^^\k) is the zeroth-order approximation to ri{k). 
Approximating (10) at the lowest order, we have 



Giz):^{kl-z')- 



(30) 



and taking into account the factorization (22) of the 
Green's function and that r]{z) has a cut in the nega- 
tive real axis, we obtain 



1 



z + ko 

Thus the second-order energy shift is 



k(2) ^ _^ V 



V ^k + ko 

kj 



/X • fkj(rA) 



(31) 



(32) 



where we have used (2) and have explicitly indicated the 
dependence of /kj and fkj from the oscillator's position 
with respect to the conducting plate. This results 
coincides with (28), of course. 

Expression (28) coincides with the well-known result 
obtained by perturbation theory at the lowest significant 
order and in the continuum limit {V = E^ — > oo), it 
yields the known atom- wall interaction energy as in 
the near zone {ta <C k^^) and as in the far zone 
{ta ^ ^(7 ) [5, 7]. In a quasi-static approach, where 
the kinetic energy of nuclei is neglected and translational 



degrees of freedom are decoupled from the electronic co- 
ordinates, the atom-plate Casimir-Polder force can then 
be obtained as minus the derivative of the energy shift 
with respect to the oscillator position Fcp = —dEg/dd. 
This result can be also generalized to more complicated 
boundary conditions by using the appropriate mode func- 
tions in the field operators. 

The expression (19) of Eg may appear more compli- 
cated than (24). However, we feel that (19) could in 
general be applied more easily because it involves G{z), 
for which we have explicit exact and approximated ex- 
pressions, rather than the function r]{z). 



IV. CONCLUDING REMARKS 

We have considered a quantum harmonic oscillator in- 
teracting with the quantum electromagnetic field in the 
continuum limit and in the presence of boundary condi- 
tions. By an improper Bogoliubov-type transformation 
we have obtained a diagonal form of the Hamiltonian in 
terms of new bosonic operators with an energy shift due 
to the oscillator-field coupling. Counter-rotating terms 
arc included in the Hamiltonian. Exact expressions for 
the ground-state energy shift have been obtained us- 
ing two different methods. As a specific case, we have 
considered an infinite conducting plate and the result- 
ing Casimir-Polder potential energy between the oscilla- 
tor and the plate, obtained from the dependence of the 
ground-state energy shift from the oscillator-plate dis- 
tance. 

We wish to stress that the validity of our method, how- 
ever, is not restricted to this specific situation, but holds 
for any boundary conditions yielding a continuous spec- 
trum for the field modes. The particular boundary con- 
dition, in fact, enters in the form of the coupling con- 
stant /kj only, and the diagonalization of the Hamilto- 
nian given in Section II, as well as the expression of the 
ground-state energy Eg given in (19) and (24), are valid 
for any form of the coupling constant. From the energy 
shift of the ground state of the interacting system, an ex- 
act expression of the Casimir-Polder force acting on the 
oscillator, due to the presence of the boundary condition, 
can be obtained. In a quasi-static approach, this imme- 
diately yields the Casimir-Polder force on the oscillator. 
We have also shown that, in the case of a perfectly con- 
ducting plate and in the lowest-order approximation in 
the coupling constant, the well-known expression for the 
second-order atom-wall Casimir-Polder energy is recov- 
ered. 

The use of the nonperturbative approach presented 
in this paper could be particularly relevant in all cases 
of the interaction of a discrete system with a quantum 
field, whenever perturbation theory cannot be applied, 
not necessarily in the framework of quantum optics or 
quantum electrodynamics. Examples are systems with 
large coupling constant and/or large density of states for 
the field modes. An example is the study of radiative 
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processes of atoms or quantum dots in a photonic crys- 
tal, when relevant transition frequencies are close to the 
band edge of the crystal. It is in fact known that in 
this case the density of states increases and has a di- 
vergence at the band-edge, resulting in the necessity of 
a nonperturbative approach. In the framework of con- 
densed matter physics, a similar situation occurs in the 
case of impurities in a quantum wire, as a consequence of 
the van Hove singularity of the electron density of states. 
We shall discuss applications to these systems in a future 
publication. 
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